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SUMMARY 


A  study  has  been  made  of  nitric  oxide  formation  in  a  laminar  CO-air  diffu¬ 
sion  flame  over  a  pressure  range  from  1  to  50  atm.  The  carbon  monoxide  (CO) 
issued  from  a  3. 06-mm-diameter  port  coaxially  into  a  coflowing  stream  of  air 
confined  within  a  20. 5-mm-diameter  chimney.  Nitric  oxide  concentrations  from 
the  flame  were  measured  at  two  carbon  monoxide  (fuel)  flow  rates:  73  standard 
cmVmin  (seem)  and  146  seem. 

As  the  pressure  was  increased  above  i  atm,  the  flames  changed  shape  from 
wide  and  convex  at  the  lowest  pressures  to  slender  and  concave  at  50  atm, 
with  the  most  noticeable  change  occurring  below  about  10  atm.  Flame  height 
decreased  as  the  pressure  was  increased  from  i  atm  to  about  1 0  atm,  but  height 
was  independent  of  pressure  at  higher  pressures.  The  regimes  of  stable  burning 
decreased  with  increasing  pressure  up  to  50  atm,  with  the  greatest  reduction 
occurring  between  1  and  20  atm. 

The  molar  emission  index  (the  moles  of  nitrogen  oxides  formed  per  mole  of 
CO  consumed)  increased  as  the  pressure  was  increased  above  1  atm,  reached  a 
maximum  at  approximately  28  to  30  atm,  and  thereafter  decreased  slowly  up  to 
50  atm. 

Overall  average  rates  of  formation  of  nitrogen  oxides  were  found  to  be 
greater  than  nitric  oxide  formation  rates  predicted  by  the  conventional  two- 
step  Zeldovich  mechanism  when  oxygen  atoms  are  assumed  to  be  in  equilibrium 
with  oxygen  molecules,  but  less  than  those  rates  predicted  by  assuming  oxygen 
atoms  to  be  in  dynamic  steady  state  with  the  CO-air  reactions.  Circumstantial 
evidence  is  offered  to  support  the  proposition  that  the  concentrations  of 
nitrogen  oxides  in  the  flames  are  at  their  equilibrium  values  for  pressures 
above  about  20  atm.  The  case  for  equilibrium  is  not  as  clear  at  the  lower 
pressures.  It  is  suspected  that  reaction  kinetics  probably  influence  or  pos¬ 
sibly  control  the  formation  of  nitrogen  oxides  at  the  lower  pressures. 

Comparison  of  the  present  data  with  data  in  the  literature  for  a  methane- 
air  diffusion  flame  shows  that  for  flames  of  comparable  flame  height  (8  to 
10  mm)  and  pseudoequivalence  ratio  (0.162),  the  molar  emission  index  of  a  CO- 
air  flame  is  significantly  greater  than  that  of  a  methane-air  flame. 


INTRODUCTION 

When  fuels  are  burned  with  air,  nitrogen  oxides  are  typically  formed  from 
the  interaction  of  nitrogen  and  oxygen  species  at  the  high  temperatures  associ¬ 
ated  with  the  combustion  process.  Some  combustion  devices,  such  as  aircraft 
turbine  engines  and  internal  combustion  engines,  operate  over  a  range  of  pres¬ 
sures  rather  than  at  one  single  pressure.  Yet,  the  effect  of  pressure  on  the 
formation  of  these  nitrogen  oxides  NOg  is  not  well  understood.  Hence,  there  is 
a  need  to  determine  experimentally  the  possible  effects  of  combustion  pressure 
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on  these  emissions.  Such  information  is  of  interest  not  only  because  of 
the  desire  to  know  the  direction  and  magnitude  of  the  effects  which  can  be 
expected,  but  also  because  such  information  can  be  used  as  a  data  base  for 
the  evaluation  of  combustion  models. 

An  effect  of  pressure  on  the  emissions  of  NOy  from  a  laminar  methane- 
air  flame  has  already  been  demonstrated  in  reference  1.  For  this  flame,  an 
increase  in  pressure  above  1  atm  is  reported  to  cause  first  an  increase  and 
then  a  decrease  in  the  molar  emission  index  (the  moles  of  NOy  produced  per  mole 
of  fuel  consumed).  The  range  of  pressures  in  that  investigation  was  from  1  to 
SO  atm,  with  the  peak  emission  occurring  at  about  <1  atm.  A  complicating  uncer¬ 
tainty  with  these  results,  however,  is  that  increasingly  large  concentrat ions 
of  carbon  particles  formed  in  the  flame  as  the  pressure  was  increased.  These 
carbon  particles  may  perturb  the  NOy-forming  reactions  from  those  which  domi¬ 
nate  when  particulate  carbon  is  not  present.  As  one  means  of  avoiding  this 
complication,  fuels  which  do  not  form  carbon,  such  as  carbon  monoxide  (CO)  and 
hydrogen  ( H 2 ) >  could  be  investigated  over  a  similar  range  of  pressures  to 
determine  whether  the  observed  maximum  in  the  molar  emission  index  is,  in  fact, 
related  to  the  presence  of  carbon  particles. 

The  present  report  details  the  results  obtained  from  an  experimental  mea¬ 
surement  of  the  molar  emission  index  from  a  CO-air  diffusion  flam*'  over  the 
pressure  range  from  1  to  SO  atm.  As  in  reference  l,  the  diffusion  flame  was 
selected  for  initial  study  of  the  effect  of  pressure  on  NOy  formation  both 
because  it  is  representative  of  the  type  of  combustion  encountered  in  many 
practical  combustion  devices  (aircraft  turbine  engines,  internal  combustion 
engines)  and  because  it  is  much  easier  to  study  safely  at  high  pressures  than 
are  premixed  flames  (no  tendency  to  flashback!.  This  investigation  was  con¬ 
ducted  in  the  same  high-pressure  facility  as  was  used  for  the  investigation  of 
reference  1.  A  qualitative  explanation  of  the  observed  results  is  given. 


SYMBOLS 

A  cross-sectional  area  for  flow  within  quartz  shield,  cm-1 

A^  proiected  area  of  base  of  flame,  cm7 

Af  surface  area  of  flame,  cm7 

concentration  of  species  i,  molos/cm, 

Cp  molar  heat  capacity,  cal/mole-K 

d  thermoco  pie  bead  diameter,  cm 

h  convection  heat-transfer  coefficient,  cal /cm7-sec-K 

T  emission  index,  (g  of  NOy  as  if  all  were  NO?) /(kg  of  CO  consumed) 

Im  molar  emission  index,  (moles  of  NOy  produced) / (mole  of  CO  consumed) 


Ira> eq  molar  emission  index  for  N0X  at  equilibrium, 

(moles  of  NOx  produced) /(mole  of  CO  consumed) 


*0 

ki 

L 

™i 

n 

• 

n 

P 

R 

rD 

rG 

ri 

rN0 

-E 

rC0 

-E 

rN0x 

rT 


equii ibrium  constant  for  oxygen  dissociation  reaction 


reaction  rate  constant  for  ith  reaction,  cm3/mole-sec  for  two-body 
reactions,  cm6/mole2-sec2  for  three-body  reactions  (number  in 
subscript  refers  to  equation  number) 


radiation  path  length,  mm 
mass  flow  rate  of  species  i,  g/sec 
total  mass  flow  rate  in  quartz  shield,  g/sec 
moles  of  species  i,  moles 
molar  flow  rate  of  species  i,  moles/sec 
pressure,  atm  (1  atm  =  101.3  kPa) 
universal  gas  constant,  82.06  cm3-atm/mole-K 
diffusional  rate  of  CO  into  flame  zone,  moles/cm3-sec' 


SiTf" 

W0  ^-ced 

Unai^0'  lon 


theoretical  estimate  of  CO  consumption  rate  based  on  global  rate 
expression  (see  appendix  D) ,  moles/cm3-sec 


reaction  rate  for  formation  of  species  i,  moles/cm3-sec 
rate  of  formation  of  NO  by  Zeldovich  mechanism,  moles/cm3-sec 


experimental  average  rate  of  consumption  of  CO,  moles/cm3-sec 
experimental  average  rate  of  formation  of  N0X,  moles/cm3-sec 


theoretical  estimate  of  CO  consumption  rate  based  on  elementary 
reaction  steps  (see  appendix  D) ,  moles/cm3-sec 


T  temperature,  K 

Ta<j  adiabatic  flame  temperature,  K 

Tb  burner  temperature  at  exit,  K 

Tc  thermocouple  temperature,  K 

Tf  flame  temperature,  K 

Tjm  temperature  required  to  force  agreement  between  Im, eq  an<1  Tm'  K 

Tref 


reference  temperature,  K 


wall  temperature  within  chamber,  K 
reaction  volume,  cm3 

volumetric  flow  rate  at  standard  temperature  and  pressure,  293  K 
and  1  atm,  standard  cm3/min 

velocity,  cm/sec 


stoich 

xo.eq 


mole  fraction  of  species  i 

equilibrium  oxygen-atom  mole  fraction  for  stoichiometric  premixed 
combustion 


*0,req  oxygen-atom  mole  fraction  required  by  equation  (4)  to  give  the 
observed  N0X  formation  rate 

AHrx  heat  of  combustion  of  CO,  cal/g 

f  emissivity 

X  thermal  conductivity,  cal/cm-K-sec 

U  viscosity,  g/cm-sec 

p  density,  g/cm3 

a  Stefan-Boltzmann  constant,  1.355  *  10-1^  cal/cm^-K^-sec 

T  residence  time,  msec 

ifc  equivalence  ratio 

4>'  pseudoequivalence  ratio 

Subscripts: 


actual  conditions 
combustion  product  effluent 
flame  conditions 
ith  species  or  ith  reaction 
flowing  into  flame  zone 

based  on  assumption  of  oxygen  atoms  in  equilibrium 
with  molecular  oxygen 

based  on  assumption  of  nonequilibrium  steady-state  balance 
for  oxygen  atoms  in  flame 
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flowing  out  of  flame  zone 


stoich  stoichiometric  conditions 


reverse  of  reaction  as  written 


Chemical  symbols 


methane 


carbon  dioxide 


hydrogen  atom 


hydrogen  molecule 


arbitrary  third  body 


nitric  oxide 


nitrogen  dioxide 


total  nitrogen  oxides  (NO  and  NO2) 


oxygen  atom 


oxygen  molecule 


hydroxyl  radical 


Abbreviations 


hydrocarbon 


parts  per  million 
standard  cmVmin 


inside  diameter 


o.d . 


outside  diameter 


dtam  diameter 

A  bar  over  a  symbol  denotes  time  averaqe. 


DESCRIPTION  OF  SYSTFM 

The  high-pressure  flame  system  used  in  the  present  study  is  the  same  as 
that  described  in  detail  in  reference  1.  It  consists  of  a  diffusion  flame 
burner  mounted  within  a  pressure  chamber,  with  provisions  for  withdrawing  com¬ 
bustion  products  and  for  viewing  and  photographi ng  the  flame.  The  burner  is  a 
3 .06-mm-diameter  tube  from  which  the  fuel,  CO,  issues  coaxially  into  a  coflow¬ 
ing  stream  of  air,  which  in  turn  is  contained  within  a  20 . S-mm-di ameter  cylin¬ 
drical  quartz  tube  concentric  with  the  burner  tube.  (See  fig.  l .)  With  a 
greater  than  stoichiometric  flow  rate  of  air,  it  is  possible  to  establish  on 
this  burner  arrangement  a  confined  laminar  diffusion  flame,  referred  to  as  a 
Burke  and  Schumann  flame  (ref.  21.  ’’’his  is  the  tvpical  over  vent  i  1  ated  diffu¬ 
sion  flame,  similar  in  shape  to  a  candle  flame. 

To  designate  the  overall  fuel-to-air  ratio,  a  pseudoequi valence  ratio  40 
is  defined  (in  keeping  with  ref.  M  as  the  actual  molar  ratio  of  fuel  to  air 
divided  bv  the  stoichiometr ic  molar  ratio  of  fuel  to  air.  t^o  term  "pseudo¬ 
equivalence  ratio"  (rather  than  "equivalence  ratio"  only)  is  used  to  draw 
attention  to  the  fact  that  the  concept  of  equivalence  ratio  does  not  have  the 
same  significance  for  a  diffusion  flame  as  it  does  for  a  premixed  flame.  In  an 
over vent il ated  diffusion  flame  (i.e.,  4*'  <11,  it  is  generally  assumed  that 

the  effective  average  equivalence  ratio  in  the  flame  zone  where  reaction  occurs 
is,  in  reality,  verv  close  to  unity. 

A  schematic  diaqram  of  the  hiqh-pressure  chamber,  burner,  and  sample  col¬ 
lection  system  is  shown  in  figure  2.  Tqnition  is  accomplished  bv  means  of  a 
resi stance-heated  hot-wire  iqniter  which  is  retractable  throuqh  a  chevron  pack¬ 
ing  gland  at  the  top  of  the  chamber.  To  iqnite  the  flame,  the  iqniter  is 
lowered  to  a  position  ’0  to  20  mm  above  the  burner  and  electrical  power  is 
applied.  Following  ignition,  the  igniter  is  retracted  out  of  the  flow  of 
combustion  products.  The  combustion  products  exit  from  the  top  of  the  quartz 
tube  (or  chimney)  and  are  swept  upward  into  a  total -sample  collector  by  the 
nitrogen  gas  used  to  pressurize  the  svstem. 

The  diluted  combustion  products  flow  out  of  the  system  throuqh  a  *>-lim 
particle  filter  and  through  a  hot-air-heated  back-pressure  regulator  used  to 
maintain  the  chamber  pressure  at  the  desired  level.  At  the  exit  of  the  Kick- 
pressure  regulator,  the  sample  enters  a  sample  transfer  line  at  near  atmo¬ 
spheric  pressure  and  flows  to  a  gas  analysis  console.  All  surfaces  contacted 
by  the  sample,  from  the  chimney  exit  to  the  qas  analysis  console,  are  of  stain¬ 
less  steel.  Color  photographs  of  the  flame  are  taken  throuqh  safety-sight 
glass  windows  located  on  opposite  sides  of  the  hiqh-pressure  chamber .  Flow 
rates  of  <30,  air,  and  Nj  are  measured  with  calibrated  linear  thermal  mass  flow¬ 
meters  having  an  accuracy  of  1  percent  of  full  scale.  Meter  ranqes  are  0  to 
200  seem  for  00,  0  to  3000  seem  for  air,  and  0  to  ^0  000  seem  for  Ni. 
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The  high-pressure  flame  system  is  operated  by  setting  the  desired  chamber 
pressure  with  the  back-pressure  regulator  and  pressurizing  with  N2.  After 
reaching  pressure,  flow  rates  for  CO,  air,  and  N2  are  adjusted  to  their  target 
values:  73  seem  or  146  seem  for  CO,  the  desired  flow  rate  for  air,  and  about 
6000  seem  for  N2.  The  flame  is  then  ignited.  Two  different  cylinders  of  CO 
were  used  in  the  experiments:  one  contained  99.3  ppm  H2  with  327  ppm  toral 
hydrocarbon  (HC)  contamination;  the  second  contained  less  than  10  ppm  H2  with 
less  than  4  ppm  HC  contamination.  The  trace  amounts  of  H2  were  intentional; 
their  purpose  was  to  serve  as  a  chain  carrier  in  initiating  and  sustaining  the 
combustion  of  the  CO.  No  difference  was  observed  in  the  data  obtained  from 
either  of  these  two  CO  cylinders.  The  air  was  zero  grade,  99.995  percent  pure 
with  an  HC  concentration  of  less  than  1  ppm;  the  N2  was  also  99.995  percent 
pure.  To  prevent  both  the  appearance  of  a  reddish-yellow  luminosity  near  the 
tip  of  the  flame  and  a  reddish-brown  deposit  on  the  burner  exit  rim  (apparently 
caused  by  the  presence  of  trace  amounts  of  iron  carbonyl  in  the  CO) ,  the  CO 
from  the  supply  cylinder  was  first  passed  through  a  dry-ice  trap  and  then 
through  an  adsorbent  bed  of  Molecular  Sieve  5A. 


EXPERIMENTAL  MEASUREMENTS 

Experimental  data  consist  of  the  chemical  composition  of  the  combustion 
products,  color  photographs  of  the  flame,  and  flame  temperatures. 


Chemical  Analyses 

Chemical  analyses  for  CO2,  CO,  O2,  and  NOx  were  performed  on  the  total  gas 
mixture  exiting  the  high-pressure  chamber.  Analyses  for  CO  and  CO2  were  per¬ 
formed  on  separate  (commercial)  nondispersive  infrared  analyzers,  and  the  anal¬ 
ysis  for  O2  was  performed  on  a  (commercial)  paramagnetic  oxygen  analyzer.  For 
the  measurement  of  NOy,  a  slightly  modified,  commercially  available  chemilumi¬ 
nescent  NO-NOx  analyzer  was  used. 

Calibration  gases  for  the  infrared  and  paramagnetic  analyzers  were  gas 
mixtures  of  CO,  CO2,  and  O2  in  N2,  with  a  stated  analytical  accuracy  of  better 
than  12  percent.  The  calibration  gases  for  the  NOx  analyzer  were  gas  mixtures 
of  NO  in  N2,  with  a  stated  accuracy  of  ±1  percent.  During  an  experimental  run, 
readings  of  the  calibration  and  zero  gases  were  frequently  interspersed  between 
readings  of  the  sample  gas  to  eliminate  possible  errors  caused  by  instrument 
drift. 


Photographs 

Color  photographs  of  the  flame  were  taken  with  a  35-mm  single-lens-reflex 
camera  fitted  wi^h  a  200-mm/f2.8  reflective  lens  and  with  color  negative  film. 
Since  the  intensity  of  the  bluish  color  of  the  flame  was  observed  to  change 
very  little  with  pressure,  a  constant  camera  exposure  of  i  /i  5  sec  was  used  at 
all  flame  conditions. 
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Flame  Temperature 


Because  of  the  difficulty  of  gaining  physical  access  to  the  flame  and 
the  steep  thermal  gradients  in  the  flame,  probing  with  thermocouples  presents 
severe  experimental  difficulties  and  was  not  attempted.  Optical  techniques 
relying  on  incandescent  carbon  particles  in  the  flame,  like  the  one  employed 
with  the  methane-air  diffusion  flame  in  reference  1,  could  not  be  used  because 
of  the  absence  of  such  particles.  Sodium  line  reversal  was  ruled  out  because 
introduction  of  sodium  into  the  flame  presents  certain  experimental  difficul¬ 
ties:  the  presence  of  sodium  in  the  system  could  perturb  the  NOy  formation 
rates  and  the  deposition  of  sodium  throughout  the  chamber  and  sample  lines  was 
to  be  avoided.  These  practical  considerations  led  to  the  decision  to  settle 
for  estimating  the  flame  temperature  at  all  pressures  from  a  knowledge  of  the 
flame  temperatures  measured  at  i  atm.  To  this  end,  the  burner  was  removed  from 
the  high-pressure  chamber  and  was  set  up  on  a  laboratory  bench  at  ambient  pres¬ 
sure.  Flames  were  established  on  the  burner  at  the  same  flow  rates  of  fuel  and 
air  as  in  the  high-pressure  experiments,  and  were  probed  with  a  chromel-alumel 
thermocouple  with  a  wire  diameter  of  0.076  mm  and  a  bead  diameter  of  0.269  mm. 
No  attempts  were  made  to  obtain  temperature  profiles  throughout  the  flames,  and 
only  peak  temperatures  were  recorded. 


RESULTS 

Experimental  measurements  were  made  at  two  flow  rates  of  CO  over  the  pres¬ 
sure  range  from  '  to  SO  atm,  73  seem  and  '46  seem.  The  limits  of  stable  burn¬ 
ing  and  the  shapes  of  the  flames  changed  markedly  with  increasing  pressure.  In 
this  section,  the  results  of  the  measurements  of  stability  and  extinction  lim¬ 
its,  flame  shape,  NOy,  and  flame  temperature  as  functions  of  pressure  are  pre¬ 
sented.  The  NOy  data  are  treated  in  the  section  "Discussion." 


Stability  and  Extinction  Limits 

A  rough  determination  of  the  stability  and  extinction  limits  was  made  with 
the  burner  prior  to  selecting  the  final  experimental  CO  and  air  flow  rate  set¬ 
tings.  These  stability  and  extinction  limits  were  determined  at  '  atm,  20  atm, 
and  SO  atm;  the  results  are  shown  in  figure  3.  Much  of  the  scatter  in  the  data 
was  caused  by  an  air  line  leak  which  was  discovered  after  the  data  had  been 
taken;  the  data  were  then  corrected  to  account  for  this  leak.  Hence,  the  indi¬ 
cated  limits  are  only  approximate.  However,  the  important  feature  to  observe 
in  figure  3  is  that  the  fuel-air  flow  rate  regimes  in  which  combustion  occurs 
decrease  considerably  as  pressure  increases.  This  phenomenon  is  similar  to 
that  reported  in  reference  '  for  a  methane-air  diffusion  flame  on  the  same 
burner  system. 

On  the  basis  of  the  data  in  figure  3,  two  fuel  flow  rates  were  selected 
for  study:  73-sccm  and  '  46-sccm  CO.  Corresponding  to  these  fuel  flow  rates, 
flow  rates  of  air  were  selected  at  1730  seem  and  2140  seem,  respectively,  giv¬ 
ing  pseudoequivalence  ratios  <!*'  of  0.100  and  0.162,  respectively.  At  73-sccm 
CO,  the  air  flow  rate  of  1730  seem  is  well  within  the  stable  region  while  at 
i  46-sccm  CO,  the  air  flow  rate  of  2' 40  seem  is  near  the  extinction  boundary  at 
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the  higher  pressures.  A  pseudoequivalence  ratio  of  0.100  could  not  be 

easily  achieved  at  the  higher  fuel  flow  rate  of  1 46-sccm  CO  because  an  air 
flow  rate  of  3476  seem  would  be  required,  a  rate  above  the  limit  of  the  air 
flow  meter. 


Flame  Shape 

Marked  changes  in  shape  and  structure  of  the  CO-air  diffusion  flames  occur 
as  the  pressure  is  increased  from  l  to  50  atm.  These  changes  can  be  easily 
seen  in  the  color  photographs  in  figure  4.  The  bottom  edge  of  each  photograph 
is  coincident  with  the  top  edge  of  the  burner,  and  the  lateral  position  of  the 
burner  is  shown  below  each  photograph.  The  bluish  color  of  the  flames  is  prob¬ 
ably  caused  by  O2  Schumann-Runge  radiation  and  CO  +  O  continuum  radiation. 
(See  ref.  3.) 

From  l  atm  to  about  5  atm,  the  flames  bow  outward  from  the  inner  fuel  core 
into  the  surrounding  air  annulus,  with  the  result  that  the  diameter  of  the 
flames  is  larger  than  the  diameter  of  the  burner  exit.  At  about  1 0  atm,  the 
diameter  of  the  flames  is  reduced,  and  the  shapes  of  the  flames  closely  approx¬ 
imate  cones.  At  pressures  above  ’0  atm,  there  is  an  increasing  concavity  to 
the  flames,  and  their  diameters  decrease  noticeably.  A  similar  behavior  was 
observed  for  the  methane-air  diffusion  flame  in  reference  1 ,  wherein  it  was 
conlectured  that  the  change  in  shape  might  result  from  the  presence  of  carbon 
particles  in  the  flame.  In  view  of  the  present  results,  this  coniecture  is 
obviously  incorrect,  and  seme  other  explanation  must  be  found.  (It  should  be 
noted  that,  according  to  classical  Burke  and  Schumann  flame  theory,  no  change 
in  the  size  or  shape  of  a  confined  laminar  diffusion  flame  should  occur  as  the 
pressure  is  changed.)  Comparison  of  the  shapes  of  the  present  flames  with 
those  of  the  methane-air  flames  discussed  in  reference  1  suggests  that  the  onlv 
effect  which  can  be  attributed  to  the  presence  of  carbon  particles  is  a  slight 
bulge  occurring  in  the  upper  third  of  the  methane-air  flames  somewhat  below 
their  tips. 

The  heights  of  the  CO-air  flames  decrease  with  pressure  from  1  atm  to 
about  tO  atm  and  thereafter  remain  relatively  constant  up  to  50  atm.  These 
flame  heights  are  plotted  in  figure  5.  The  decrease  in  the  height  of  the  flame 
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at  a  1 46-sccm  flow  rate  is  from  about  9—  mm  to  about  7—  mm,  or  roughly  20  per- 
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cent;  the  decrease  for  the  flame  at  a  73-sccm  flow  rate  is  from  about  6  mm  to 
about  3  mm,  or  roughly  50  percent. 


Oxides  of  Nitrogen 

Results  of  the  experimental  measurements  on  the  total  oxides  of  nitrogen 
are  shown  in  figure  6  as  Ip,,  molar  emission  index  of  NOy,  as  a  function  of 
pressure.  Experimentally,  this  index  is  determined  as  the  ratio  of  the  mea¬ 
sured  concentrations  of  NOy  to  CO2  in  the  effluent  from  the  high-pressure  cham¬ 
ber.  Ttie  molar  emission  index  is  analogous  to  the  emission  index  customarily 
used  to  report  data  on  gas  turbine  engines;  the  latter  is  defined  as  the  number 
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of  grams  of  NOy  produced,  as  if  all  nitrogen  oxides  were  present  as  NO2,  per 
kilogram  of  fuel  burned.  According  to  these  definitions,  the  emission  index 
for  a  CO  flame  ib  1045  times  the  molar  emission  index. 

To  obtain  the  data  points  in  figure  6,  the  CO  flow  rates  were  set  at  their 
target  values,  three  or  more  runs  were  made  at  different  air  flow  rates,  and 
concentrations  of  NOy,  CO,  CO2,  and  O2  were  measured.  Then  values  of  cal¬ 

culated  from  these  data  were  plotted  against  pseudoequivalence  ratio  1  and 
were  read  back  at  the  desired  values  of  <t> '  (appendix  A  gives  the  method  of 
calculating  $').  This  procedure  effectively  eliminated  the  need  for  setting 
the  air  flow  rates  at  precise  values,  a  difficult  task  at  times.  It  also 
avoids  a  reliance  on  the  accuracy  of  the  air  flow  meter.  On  the  whole,  over 
the  usual  ranges  of  4'*  employed  (roughly  0.08  to  O.M  for  73-sccm  CO  and  0.14 
to  0.18  for  1 46-sccm  CO),  there  was  no  pronounced  or  consistent  variation  of 
molar  emission  index  with  d>*  which  could  be  attributed  to  other  than  experi¬ 
mental  error.  Above  about  30  atm,  however,  there  did  appear  to  be  a  general 
decrease  of  Im  with  <V,  although  this  decrease  was  not  large.  For  purposes 
of  comparison,  the  molar  emission  index  for  the  methane-air  flame  studied  in 
reference  l  has  been  included  in  figure  6. 

Figure  8  shows  that  for  both  CO  flames,  the  molar  emission  index  starts 
at  a  relatively  low  value  at  low  pressure  and  increases  as  the  pressure  is 
increased.  The  molar  emission  index  for  the  flame  with  the  lower  flow  rate 
reaches  a  maximum  of  about  2.5  *  10“^,  while  the  Im  for  the  flame  with  higher 
flow  rate  reaches  a  maximum  of  about  4.6  x  1 0"-*,  both  maxima  occurring  at 
approximately  28  to  30  atm.  At  still  higher  pressures,  the  Im  for  both 
flames  levels  off  and  drops  slightly.  These  results  are  in  contrast  with  those 
for  the  CH4  flames,  which  exhibit  a  much  smaller  increase  with  pressure,  attain 
an  earlier  maximim  (at  about  9  atm),  and  then  decrease  to  near  the  low-pressure 
value  again  at  50  atm.  This  comparison  between  the  CH4  flame,  with  its  low 
flow  rate  of  42  scan,  and  the  CO  flames  with  a  higher  flew  rate  may  at  first  not 
seem  proper;  however,  it  is  noted  that  in  terms  of  flame  height  (8  to  ’0  mm), 
pseudoequivalence  ratio  (0.162),  and  air  flow  rate  (2450  seem),  the  CH4  flame 
is  quite  similar  to  the  1  46-sccm  CO  flame.  Yet  the  Im  for  this  CO  flame  is 
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roughly  4—  times  the  Im  for  the  CH4  flame.  One  obvious,  at  least  partial, 

explanation  is  that  the  CO  flames  probably  attain  higher  flame  temperatures  by 
virtue  of  their  higher  adiabatic  flame  temperatures  -  2384  K  for  CO  as  compared 
with  2230  K  for  CH4.1  Furthermore,  because  radiation  from  the  CH4  flame  is 
significant,  its  actual  flame  temperature  is  likely  to  be  even  lower  than  the 
154  K  difference  in  adiabatic  flame  temperatures  would  indicate.  At  the  lower 
pressures,  where  carbon  radiation  from  the  CH4  flames  is  not  yet  substantial, 
the  Im  values  for  both  flames  are  in  closer  agreement.  The  effects  of  tem¬ 
perature  on  Im  are  covered  in  more  detail  in  the  section  "Discussion." 


’At  1  atm.  Although  there  is  a  slight  increase  in  adiabatic  flame  temper¬ 
ature  with  increase  in  pressure  for  both  fuels,  the  increase  is  not  dissimilar. 
The  values  of  adiabatic  flame  temperature  were  calculated  using  the  chemical 
equilibrium  program  of  Svehla  and  McBride  (ref.  4). 
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Only  two  studies  could  be  found  in  the  literature  where  measurements  of 
N0X  in  CO-air  flames  at  elevated  pressures  are  reported.  These  data  are  com¬ 
pared  with  the  present  data  in  figure  7  on  a  log-log  plot.  The  data  points  A, 
B,  and  C  are  from  Haber  and  Coates  (ref.  5)  for  a  diffusion  flame  burner  simi¬ 
lar  in  design  to  the  one  used  in  the  present  study,  but  with  a  i . 5-mm-diameter 
fuel  port.  Most  of  the  Haber  and  Coates  data  are  not  for  air  but  for  support¬ 
ing  atmospheres  consisting  of  O2/N2  mixtures  of  varying  composition.  Points  A, 
B,  and  C  are  the  only  data  points  at  which  direct  comparisons  with  the  present 
results  are  possible.  Even  though  these  data  are  for  different  flow  rates  of 
CO,  and  were  obtained  with  a  burner  of  smaller  diameter,  they  are  within  rea¬ 
sonable  proximity  of  the  present  data. 

The  data  of  Newitt  and  Lamont  (ref.  6),  shown  as  point  D,  are  for  an 
opposed- jet  diffusion  flame  and  represent  their  only  data  for  air.  Although 
they  did  not  state  the  flow  rates  of  their  reactants,  it  is  reasonable  to 
assume  that  because  of  the  opposed-jet  flow,  there  must  have  been  a  turbulent, 
fairly  well-mixed  reaction  zone.  This  type  of  reaction  system  might  be 
expected  to  approach  more  closely  that  of  premixed  combustion  rather  than 
diffusional  combustion,  or,  at  the  very  least,  to  be  somewhere  between  the  two. 
However,  there  are  other  problems  with  the  data  of  Newitt  and  Lamont  that  make 
an  interpretation  of  their  data  difficult.  In  particular,  large  concentrations 
of  unburned  CO  remained  in  the  product  gases  even  at  relatively  low  values  of 
<*>',  demonstrating  that  combustion  was  far  from  complete. 


Flame  Temperature 

As  mentioned  previously,  temperature  data  were  obtained  for  flames  at  both 
flow  rates  at  i  atm  only.  With  these  data  and  estimates  of  heat  losses  from 
the  flames,  a  reasonable  trend  of  flame  temperature  with  increasing  pressure 
was  calculated.  The  raw  thermocouple  data  were  first  corrected  for  radiation 
losses  by  the  methods  detailed  in  appendix  B.  Flame  temperature  as  a  function 
of  pressure  was  then  estimated  by  the  procedure  described  in  appendix  C.  No 
attempt  was  made  to  correct  the  thermocouple  data  for  the  heat  conduction 
losses  that  can  arise  from  the  steep  temperature  gradients  in  the  flame.  (In 
severe  cases,  such  losses  could  possibly  produce  indicated  readings  as  much  as 
l 00  K  too  low,  but  these  are  extremely  difficult  to  estimate  accurately.)  The 
resulting  temperature  curves  are  shown  in  figure  8  along  with  a  plot  of  the 
adiabatic  flame  temperature  for  a  CO-air  flame  at  <t>  =  1  .0.2  These  temperature 
curves  are  believed  to  be  fairly  reasonable  estimates  of  the  true  flame  temper¬ 
atures.  For  the  t 46-sccm  CO  flame,  temperature  increases  approximately  370  K 
between  l  atm  and  50  atm;  for  the  73-sccm  CO  flame,  there  is  a  corresponding 
increase  of  approximately  480  K.  These  increases  compare  with  an  increase  of 
1 40  K  in  adiabatic  flame  temperature  over  the  same  pressure  range. 


2It  may  bear  repeating  here  that  even  though  the  overall,  or  pseudoequiva¬ 
lence  ratio  <*>'  may  be  considerably  smaller  than  unity,  the  effective  average 
equivalence  ratio  within  the  reaction  zone  must  be  quite  close  to  unity.  The 
adiabatic  flame  temperature  was  calculated  from  the  equilibrium  computer  pro¬ 
gram  of  Svehla  and  McBride  (ref.  4). 


DISCUSSION 


In  this  section,  the  averaqe  rate  of  formation  of  N0X  is  estimated  from 
the  molar  emission  index  and  the  averaqe  rate  of  consumption  of  CO.  The  ques¬ 
tion  of  the  combustion  reactions  approach inq  equilibrium  is  treated,  and  the 
averaqe  rates  of  formation  of  NOjj  are  discussed  and  compared  with  various  theo¬ 
retical  predictions.  Finally,  the  shapes  of  the  molar  emission  index  curves 
are  discussed  from  the  point  of  view  of  their  be i nq  influenced  by  approach  to 
equilibrium,  particularly  at  the  hiqher  pressures. 


Averaqe  Rate  of  Formation  of  N0X 

An  estimate  of  the  averaqe  rate  of  formation  of  N0X  within  the  riame  was 
obtained  as  follows.  At  steady-state  conditions  and  after  complete  combustion 
of  the  f uel , 
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where  npp  is  the  molar  flow  rate  of  CO  into  the  flame  and 


nNOv 


is  the 


molar  flow  rate  of  NOx  out  of  the  flame.  Aftei  tear r anqement  of  this  expres¬ 
sion  and  division  by  a  reaction  volume  VRX,  there  results 


nNOx  nOO 

vRx  m  vRx 

If  the  reaction  volume  in  which  Nf^  is  formed  is  taken  to  lv>  approximately 
equal  to  (thouqh  not  necessarily  coincident  with)  the  reaction  volume  in  which 
CO  is  consumed,  this  expression  can  I  to  rewritten  as 
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rNOx  *  Tmrai  (1) 

where  rNOx  * M  experimental  averaqe  rate  of  formation  of  NOx  and  r(X1  is 

the  experimental  averaqe  rate  of  consumption  of  CO.  An  estimate  of  the  maqni- 
-F 

tudo  of  rpn  was  obtained  tty  dividing  the  known  molar  flow  rate  of  CO  n,^  by 
the  volume  of  the  flame  zone,  oxper i mental  1 y  obtainable  f trim  the  color  photo¬ 
graphs  in  figure  4.  (These  inaction  volumes  should  lie  a  fairly  accurate  mea¬ 
sure  of  the  volumes  in  which  fuel  is  actually  consumed  since  the  blue  zones  in 
fiqure  4  are  largely  the  result  of  chem i 1 umi nescence  from  the  CO  +  0  reac¬ 
tion,  and  since  combustion  was  observed  to  !v>  virtually  complete  within  the 
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flame.)  Calculation  of  'NOs  from  equation  (1)  implies  that  the  reaction 


vo  1  ume  for  NOx  formation  is  roughly  equal  to  that  for  CO  combustion.  Although 
this  is  probably  not  strictly  true,  it  is  certainly  reasonable  to  assume  that 
these  two  volumes  are  roughly  proportional.  Since  these  reaction  volumes  prob¬ 
ably  do  not  differ  by  much  more  than  a  factor  of  2  either  way,  in  the  absence 
of  Information  to  the  contrary,  they  have  been  taken,  for  convenience,  to  tie 
equa 1 . 

The  photographs  in  figure  4  were  used  to  determine  the  magnitudes  of  the 
CO  reaction  volumes  VRx  by  tracing  the  thin  blue  reaction  /ones  at  the  sur¬ 
faces  of  the  flames  and  by  analytically  revolving  the  areas  so  defined  about 
the  center  line  of  the  flames.  The  resulting  reaction  volumes  are  shown  in 
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figures  9(a)  and  9(b).  Values  for  rpp  based  on  these  reaction  volumes 
are  plotted  as  a  function  of  pressure  in  figure  10(a).  Also  plotted  in  fig¬ 
ure  10(a)  are  two  different  predictions  of  CO  combustion  rates  r^  and  rT 
for  a  global  rate  expression  (ref.  7)  and  a  set  of  elementary  reaction  steps 
for  CO  combustion,  respectively.  Both  are  based  on  the  assumption  of  an  effec¬ 
tive  average  equivalence  ratio  of  unity  in  the  flame  /one  and  on  initial  com¬ 
position.  Details  of  the  methods  used  for  evaluating  r(;  and  rT  are  given 
in  appendix  D. 

Because  rj;  and  r^  are  based  on  initial  composition  and  an  average  <f< 
of  unity,  they  very  likely  represent  the  approximate  maximum  rates  of  combus¬ 
tion  of  CO  in  air  for  the  present  flame  conditions.  Departure  of  the  expeti- 
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mental  rates  r^  from  these  predicted  rates  could  mean  that  the  assumption  of 
an  average  ^  of  unity  is  not  a  good  one,  that  the  choice  of  initial  composi¬ 
tion  for  the  reactants  is  not  appropriate,  that  the  use  of  Tr  as  an  average 
temperature  throughout  the  flame  /one  is  not  proper,  that  diffusion  is  reallv 
controlling  the  rate  of  reaction,  or  some  combination  of  these.  iv>  treat  the 
reaction  /one  in  a  diffusion  flame  in  such  simplified  terms  as  these  obviouslv 
oversimplifies  the  true  situation.  However,  this  simplified  picture  allows  the 
inference  from  figure  10(a)  (by  virtue  of  the  fact  that  the  predicted  rates 
are,  for  the  most  part,  greater  than  the  experimental  rates)  that  the  rates  of 
combustion  of  CO  in  the  present  flames  are  largely  diffusion  controlled,  cer¬ 
tainly  at  pressures  above  2  atm  or  so. 

An  additional  reason  to  suspect  diff us  tonal  control  of  the  present  CO-air 
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flames,  even  at  all  pressures,  is  the  proximity  of  the  observed  rates  r^ 
for  flames  at  both  flow  rates.  Because  the  flames  have  significantly  dif¬ 
ferent  flame  temperatures  they  also  would  be  expected  to  have  much  different 
reaction  rates  if  reaction  kinetics  were  actually  controlling.  However,  since 
these  rates  do  not  differ  significantly,  and  since  rates  of  diffusion  are  rela¬ 
tively  insensitive  to  temperature,  it  is  reasonable  to  infer  that  the  combus¬ 
tion  rates  in  both  flames  are  largely  diffusion  controlled.  A  stronger  argu¬ 
ment  for  diffusional  control  can  he  obtained  by  estimating  the  diffusion.!! 
rates  of  00  into  the  flame  /one  rp.  This  estimate  was  carried  out  hv  the 
method  described  in  appendix  E,  and  the  resulting  values  of  ip  are  shown  in 
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figure  10(b).  (Because  these  rp  virtually  overlay  the  r^  in  figure  10(a) 


1  1 


they  have  been  plotted  on  a  separate  figure  for  ease  in  visualization.)  The 
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almost  exact  agreement  between  rpg  and  the  estimated  rates  rp  seems  con¬ 
vincing  evidence  that  the  combustion  reactions  are,  in  fact,  diffusion  con¬ 
trolled.  The  important  implication  of  the  flames  being  diffusion  controlled 
is  that  the  combustion  reactions  approach  equilibrium. 


The  average  rates  of  formation  of  nitrogen  oxides  in  the  flame  were  esti- 
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mated  using  equation  (1)  and  figures  6  and  10(a)  for  Im  and  r^p,  respec- 
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tively.  The  resulting  rates  rN0x  aro  plotted  in  figure  11,  Also  plotted  in 

figure  11  are  two  sets  of  theoretical  rate  predictions,  each  based  on  a  differ¬ 
ent  set  of  assumptions,  but  both  employing  the  conventional  two-step  Zeldovich 
mechanism  for  NO  formation: 


0  +  N2  *  NO  +  N 


(2) 


N  +  02  NO  +  O 


(3) 


By  standard  treatments  (see,  for  instance,  ref.  I),  the  rate  of  formation  of  NO 
by  the  Zeldovich  mechanism  is  given  by 

rNO  “  2k2cOcN2  (4> 

where  k2,  the  forward  reaction  rate  constant  of  reaction  (2),  is  (refs.  1,  8, 
and  9) 


k2  -  1.36  x  10'4  e~37947'/>r 

One  set  of  the  theoretical  rate  predictions  in  figure  11  is  based  on  the 
assunption  that  oxygen  atoms  are  in  equilibrium  with  molecular  oxygen,  and  the 
other  is  based  on  the  assumption  that  there  exists  superequilibr inn  steady- 
state  concentrations  of  oxygen  in  the  flames.  To  obtain  the  theoretical  rate 
predictions  for  the  situation  of  oxygen  atoms  in  equilibrium  with  molecular 
oxygen,  the  equilibrium  constant  (refs.  1  and  9) 


Kq  ■ 


(Co)2 

~^T 


3.56  x  104  -  e-59386/*1, 
T 


is  substituted  into  equation  (4)  to  eliminate  the  unknown  oxygen-atom  concen¬ 
tration  Cq.  In  terms  of  mole  fractions,  this  rate  expression  then  becomes 


rN0  -  6.90  x  10’ 3  — -  e“67640/'r  xo£2xM2 
T2 


(1) 
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0.21 


This  expression  was  evaluated  at  the  maximum  possible  values  Xq, 


and 


Xfj2  ■  0.79,  at  T  «  Tf,  and  is  labeled  (rNO^O,eq*  Rates  so  obtained  are  the 


maximum  rates  of  formation  of  NOy  attainable  when  oxygen  atoms  are  in  equilib¬ 
rium  with  molecular  oxygen  since  the  effective  mole  fraction  for  oxygen  in  the 
zone  where  nitric  oxide  is  formed  must  be  substantially  less  than  the  0.21 
used,  and  the  effective  mole  fraction  for  nitrogen  is  probably  somewhat  less 
than  the  0.79  used.  For  instance,  in  reference  1  it  was  estimated  that  in  a 
methane-air  diffusion  flame,  a  reasonable  value  for  the  effective  oxygen  mole 
fraction  is  on  the  order  of  0.04.  If  this  were  also  the  case  for  the  present 
CO-air  diffusion  flames,  the  predicted  rates  <rNo'o,eq  *n  f*9ure  11  would  be 
about  0.4  times  those  indicated  -  lower  yet,  and  farther  away  from  agreement 


with  the  experimentally  observed  average  rates  rNOx 


Even  though  these  exper¬ 


imentally  observed  rates  are  rough  estimates,  it  seems  unlikely  that  they  could 
be  in  error  by  a  factor  as  large  as  100  or  more,  which  they  would  have  to  be  at 
1  atm  if  the  predicted  rates  (tNO^O,eq  are  correct. 


One  of  the  more  reasonable  explanations  for  the  higher  experimental  rates 
is  that  the  oxvqen  atoms  in  the  flame  zone  are  not  in  equilibrium  with  the 
molecular  oxygen  but  are  present  in  concentrations  greatly  in  excess  of  equi¬ 
librium  values.  Such  an  explanation  is  consistent  with  the  observation  that 
the  experimental  rates  approach  the  predicted  rates  (rN0*0,eq  at  higher 

pressures  where  superequ i l i hr ium  oxygen-atom  concentrations  are  less  likely. 
Calculations  of  the  magnitude  required  for  these  superequi 1 ibr i urn  oxygen-atom 
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concentrat ions  for  the  predicted  rates  to  agree  with  rNOy  can  ^  made  by  sub- 

stituting  rNOx  f°r  rNO  *n  equation  (4),  takinq  t  =  Tf  and  xN_,  =  0.79,  and 

solving  for  the  oxvgen-atom  concentration.  i’he  resulting  concentrations  are 

stoi ch 

plotted  in  figure  12  as  the  ratio  xo,req/xO,eq  >  where  the  denominator  is  the 
equilibrium  oxygen-atom  mole  fraction  for  stoichiometric  premixed  combustion  at 
Tf.  Extremely  high  values  of  the  required  oxygen-atom  ratio  occur  at  the  lower 
pressures,  but  they  decrease  continually  as  the  pressure  is  increased. 


no  obtain  a  more  realistic  estimate  of  the  oxygen-atom  concentrations,  an 
oxygen-atom  balance  was  made.  This  balance  is  based  on  the  assumption  that  a 
steady-state  condition  exists  in  the  flame  wherein  the  rate  of  generation  of 
oxygen  atoms  is  exactly  counterbalanced  by  the  rate  of  removal  of  oxygen  atoms. 
The  details  of  this  approach  are  given  in  appendix  F.  Bv  this  method  oxygen- 
atom  concentrations  have  been  calculated  and  converted  to  the  ratio 

f?  to  i  ch 

xO  ns/xO  eq  and  are  plotted  in  figure  12.  For  the  1 46-sccm  00  flame,  this 
concentration  of  oxygen  atoms  is  considerably  greater  than  that  required  to 
produce  the  experimentally  observed  rates,  and  for  the  73-sccm  00  flame,  this 
concentration  is  qreater  than  that  required  for  all  pressures  above  about 
3  atm.  For  pressures  below  about  3  atm,  the  estimated  oxygen-atom  concentra¬ 
tion  is  less  than  that  required,  being  as  much  as  an  order  of  magnitude  less  at 


1  atm.  This  lower  value  of  xo>gg  for  the  73-sccm  CO  flame  below  3  atm  implies 
that  there  is  an  insufficient  concentration  of  oxygen  atoms  generated  by  the 

-E 

combustion  process  below  3  atm  to  account  for  the  observed  rate  fNOx  bY 
Zeldovich  mechanism.  However,  because  of  the  size  of  the  potential  errors 


— E 

involved  in  calculating  r*jox»  selecting  an  appropriate  value  for  xN^,  in 

making  the  steady-state  analysis  for  xq  gs,  as  well  as  in  estimating  Tf,  there 
is  insufficient  cause,  on  these  grounds  alone,  to  invalidate  the  Zeldovich  mech¬ 
anism  as  being  the  basic  mechanism  for  the  production  of  nitric  oxide. 

Rates  of  formation  of  nitric  oxide  based  on  these  steady-state  oxygen-atom 
concentrations  Uno)o,ss  have  been  calculated  using  equation  (4)  and  are  com- 

—E 

pared  with  the  experimentally  observed  rates  rN0x  in  figure  11.  For  this  cal¬ 
culation,  the  nitrogen  concentration  was  taken  to  be  xN.,  =  0.556  (i.e.,  that 

nitrogen  concentration  existing  in  a  stoichiometric  mixture  of  CO  and  air). 

These  rates,  (tNO^O,ss*  are  not  only  considerably  greater  than  the  rates  pre¬ 
dicted  on  the  basis  of  oxygen-atom  equilibrium  but  are  also  greater  than  the 
experimentally  observed  rates  except  for  the  73-sccm  CO  flame  below  about  3  atm. 
Clearly,  the  reason  why  (rNo)o  ss  is  considerably  greater  than  (rNO^O.eq  ir' 
that  the  oxygen-atom  concentrations  based  on  the  steady-state  analysis  are  con¬ 
siderably  greater  than  those  based  on  the  equilibrium  assumption.  (Examination 
of  the  analysis  in  appendix  F  shows  that  this  is  the  case  because  the  reaction 
step  with  the  fastest  rate  of  production  of  oxygen  atoms  is  the  combustion 
reaction 


CO  +  C>2  +  CO  2  +  O 
and  not  the  dissociation  reaction 
02  +  M  +  20  +  M 

which  it  would  have  to  be  for  atomic  oxygen  to  be  in  equilibrium  with  molecular 
oxygen . ) 

From  figure  11,  it  is  apparent  that  if  the  Zeldovich  mechanism  is  the  cor¬ 
rect  mechanism  for  nitric  oxide  formation  in  CO-air  diffusion  flames,  and  if 

_E 

the  experimentally  determined  rates  rN0x  arG  approximately  correct,  then  the 

true  average  oxygen-atom  concentrations  in  the  flames  are  not  as  large  an  those 
calculated  on  the  basis  of  the  steady-state  analysis,  that  is,  on  the  basis  of 
purely  kinetic  considerations.  This  is  certainly  the  case  at  pressures  above 

-E 

about  15  atm  or  so,  where  (rN0>0,ss  sufficiently  greater  than  rj*)x  so 
that  it  is  unlikely  that  the  difference  can  be  attributed  to  experimental 
16 
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error.  Below  15  atm,  and  particularly  below  about  3  or  4  atm,  the  case  is  not 
so  clear.  Precisely  the  same  conclusion  can  be  reached  by  inspection  of  fig- 

ure  12,  where  xo,req/*0,eq  *s  less  than  Xq  sa/xo,eq  at  except  the 
lower  pressures.  This  result  leads  one  to  consider  the  possibility  that  per¬ 
haps  reaction  kinetics  may  not  be  controlling  the  rate  of  formation  of  nitric 
oxide,  but  that  diffusion  might  be  playing  a  role.  If  diffusion  were,  in  fact, 
controlling,  then  it  necessarily  follows  that  the  NOy  concentrations  in  the 
flames  should  approach  their  equilibrium  values.  This  eventuality  is  consid¬ 
ered  in  the  following  section. 


Molar  Emission  Index 

Equilibrium  values  of  the  molar  emission  index  calculated  for  stoichiomet¬ 
ric  combustion  and  t  =  Tf  are  plotted  in  figure  1 3  as  lm#(»q  »l°ng  with  the 
curves  of  Tro  from  figure  6.  The  similarity  in  shape  between  these  curves  is 
obvious.  However,  lm,eq  rises  faster  and  peaks  sooner  than  does  Im.  Above 
about  28  atm,  the  curves  of  !m,eq  Parallel  those  of  Tm  fairly  well  but  are 
roughly  20  to  25  percent  lower.  Yet,  the  general  similarity  in  shape  suggests 
the  possibility  that  the  NOy  concentrations  in  the  flames  might  be  approaching 
their  equilibrium  values,  particularly  at  the  higher  pressures  where  the  curves 
of  and  Tm^q  are  reasonably  parallel,  and  where  the  required  oxygen-atom 

ratios  xo,req/xo*eqCh  *n  fi9yre  12  are  decreasing  to  lower,  not  so  excessive 
values.  At  pressures  above  28  atm,  separate  calculations  indicate  that  an 
increase  in  Tf  of  only  about  50  K  for  both  flames  would  be  required  to  bring 
Tm,eq  into  exact  agreement  with  Tm,  a  variation  certainly  well  within  the 
error  of  estimating  Tf.  Circumstantial  evidence  of  a  similar  nature  (i.e., 
the  agreement  of  ^m,eq  an<*  Tm'  was  in  reference  1  for  a  CH4-air  dif¬ 

fusion  flame  in  support  of  the  proposition  that  the  NOy  concentrations  in  that 
flame  had  approached  equilibrium  levels  at  pressures  above  about  20  atm  or  so. 
Unlike  the  present  OO-air  flames,  however,  it  was  clear  for  the  CH4~air  flame 
that  equilibrium  levels  of  NOy  had  not  been  approached  at  the  lower  pressures, 
and  that  kinetics  were  still  controlling.  For  the  present  CO-air  flames 
the  case  is  not  so  clear.  However,  because  of  the  very  large  values  of 

stoich 

xO,req/xO,eq  at  the  lower  pressures  (see  fig.  12),  it  is  suspected  that 
kinetics  are  controlling  even  though  reasonably  close  to  Tm. 

Because  the  experimentally  estimated  temperatures  Tf  are  approximate  in 
nature  and  because  equilibrium  levels  of  nitric  oxide  are  quite  sensitive  to 
temperature,  additional  calculations  were  made  to  determine  those  temperatures 
which  would  have  had  to  exist  in  the  flames  in  order  for  Tm,eq  to  coincide 
identically  with  Tm  over  the  entire  range  of  pressures.  These  temperatures 
are  plotted  in  figure  1 4  as  TT  .  Comparison  of  these  temperatures  with  Tf 

shows  the  average  deviation  between  them  to  be  about  40  K,  a  value  well  within 

the  experimental  error  of  the  estimated  temperatures  Tf.  Such  close  agreement 

between  T»  and  Tf  is,  nevertheless,  circumstantial  and  therefore  insufft- 
*111 

cient  to  prove  that  the  NOy  concentrations  approach  their  equilibrium  values  at 
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all  pressures.  Further,  because  of  the  uncertain  nature  of  Tf  anti  the  great 
sensitivity  of  Im,eq  Tf*  case  for  or  against  equilibrium  NO*  cannot 

he  resolved  simply  by  consideration  of  the  molar  emission  index.  At  the  lower 

pressures,  the  very  large  values  of  *o,  req//*<%  eq°h  tend  to  argue  aqainst 
equilibrium  there.  Yet,  the  possibility  that  NO*  concentrations  approach  equi¬ 
librium  within  the  flames  at  the  higher  pressures  ahove  about  20  atm  seems  dis¬ 
tinctly  real. 

To  rule  out  the  possibility  that  the  shape  of  the  lm  curves  might  be 
caused  by  changes  in  residence  time,  residence  times  were  calculated  from  the 
expression 


vRx 

T  *  - 

V*„1 

where  Vrx  is  the  reaction  volume  (plotted  in  figs.  R(a)  and  <>(M1  and 
is  the  volumetric  flow  rate  of  00  and  air  at  an  average  equivalence  ratio  of 
unity  and  at  T  *  Tf.  Resulting  values  of  T  are  shown  in  figure  IS,  where  T 
first  decreases,  then  passes  through  a  minimum  between  1  and  7  atm,  and  there¬ 
after  continues  to  increase  with  increasing  pressure.  Since  this  trend  is  con¬ 
trary  to  the  one  for  Tm,  it  is  obvious  that  the  shape  of  the  Tm  curves  is 
not  determined  by  residence-time  considerations?  in  particular,  the  decrease  in 
Tm  observed  at  the  higher  pressures  is  not  caused  by  a  corresponding  decrease 
in  1  . 

To  help  resolve  this  question  of  the  NOy-forminq  reactions  approaching 
equilibrium,  a  mathematical  model  of  a  CO-air  diffusion  flame  would  be 
extremely  useful.  tf  such  a  model  employed  reasonable  kinetic  mechanisms 
describing  combustion,  and  if  the  model  adequately  predicted  the  observed 
chanqes  in  flame  structure  with  pressure,  then  there  would  lie  some  hope  that 
the  basic  details  of  the  NOy  formation  processes  could  also  lx*  adequately 
modeled,  and,  hence,  better  understood.  Recently,  a  variety  of  analytical 
approaches  and  computer  codes  capable  of  mathemat ical ly  model inq  both  the  com¬ 
plicated  fluid  dynamic  and  chemical  aspects  of  diffusion  flames  have  hoeome 
available  (refs,  'fl  to  1*11,  and  these  could  possibly  prove  useful. 

As  a  final  comment,  it  was  speculated  in  reference  i  that  among  the  vari¬ 
ous  factors  which  could  possibly  account  for  the  presence  of  the  maximum  in  tm 
for  a  O^-air  diffusion  flame  is  the  large  amount  of  carbon  in  the  flame.  The 
present  results  clearly  show  that  the  presence  of  carbon  in  the  flame  is  not  a 
requirement  for  such  a  maximum.  Although  such  carbon  could  still  conceivably 
affect  the  location  of  the  maximum,  it  is  obviously  not  necessary  for  such  a 
maximum  to  exist.  Based  on  the  present  data  and  analysis,  the  preferred 
explanation  for  the  shape  of  the  Tm  curves  as  a  function  of  pressure  for  both 
the  PO-air  and  CH4-air  diffusion  flames  is  that  equilibrium  levels  of  NOy  are 
approached  in  the  flames  at  the  higher  pressures. 
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CONCLUDING  REMARKS 


A  study  has  been  made  of  nitric  oxide  formation  in  a  laminar  CO-air  diffu¬ 
sion  flame  over  a  range  of  pressures  from  l  to  SO  atm.  The  carbon  monoxide 
(fuel)  issued  from  a  3. 06-mm-diameter  port  coaxially  into  a  coflowing  stream  of 
air  confined  within  a  20. 1-mm-diameter  chimney.  Experimental  measurements  of 
nitric  oxide  concentrations  were  made  at  two  carbon  monoxide  (fuel)  flow  rates: 
73  standard  cm^/min  (seem)  and  MS  seem.  The  flow  rates  were  held  constant 
over  the  entire  range  of  pressures  from  1  to  so  atm.  The  stability  and  extinc¬ 
tion  limits,  flame  shape,  flame  heiqht,  NCl^,  and  flame  temperature  were  all 
affected  by  changes  in  pressure.  Rased  on  these  experimental  data  and  their 
analysis,  the  following  observations  and  conclusions  are  offered: 

1.  The  shape  of  CO-air  diffusion  flames  chanqes  from  wide  and  convex  at 
1  atm  to  slender  and  concave  at  SO  atm.  The  most  noticeable  change  in  shape 
occurs  below  about  ’0  atm. 

2.  The  height  of  the  flames  decreases  with  pressure  from  l  atm  to  about 

10  atm  or  so,  and  thereafter  remains  relatively  constant  up  to  SO  atm.  For  the 

1  1 

flame  with  i 46-sccm  flow  rate,  the  decrease  is  from  about  mm  to  about  7-  mm; 

2  2 

for  the  73-sccm  flame,  the  decrease  is  from  about  6  mm  to  a  little  over  3  mm. 

3.  The  regimes  of  stable  burning  decrease  as  the  pressure  is  increased, 
with  the  greatest  reduction  occurring  between  l  atm  and  20  atm. 

4.  The  concentrat ions  of  nitrogen  oxides  produced  by  the  flames  wore  mea¬ 
sured  and  reported  as  molar  emission  index  (the  moles  of  nitroqon  oxides  formed 
per  mole  of  carbon  monoxide  consumed).  For  both  flames,  the  molar  emission 
index  increases  as  the  pressure  is  increased  above  i  atm,  reaches  a  maximum  in 
the  neighborhood  of  28  to  30  atm,  and  thereafter  decreases  slowly  up  to  *30  atm. 
Hence,  at  all  pressures  above  i  atm  and  up  to  at  least  10  atm,  the  molar  emis¬ 
sion  index  is  greater  than  that  at  i  atm. 

5.  Comparison  of  the  present  data  with  data  in  the  literature  for  a 
methane-air  diffusion  flame  shows  that  for  flames  of  comparable  flame  heiqht 
(8  to  10  mm)  and  pseudoequivalence  ratio  (0.162),  the  molar  emission  index  of 

1 

a  CO-air  flame  can  be  as  much  as  4—  times  that  of  a  methane-air  flame. 

2 

6.  Those  few  data  in  the  literature  on  nitric  oxide  concentrat ions  gen¬ 
erated  by  coaxial  CO-air  diffusion  flames  are  in  reasonable  agreement  with  the 
present  data. 

7.  Overall  average  reaction  rates  of  carbon  monoxide  were  determined  for 
the  diffusion  flames.  Comparison  of  these  rates  with  kinetic  predictions  and 
diffusion  rates  suggests  that  the  carbon  monoxide  combustion  reactions  are  dif¬ 
fusion  controlled.  The  implication  of  this  is  that  these  combustion  reactions 
approach  equilibrium. 


1  <» 


8.  Overall  average  rates  of  formation  of  nitrogen  oxides  were  estimated 
from  the  data.  These  rates  are  greater  than  nitric  oxide  formation  rates  pre¬ 
dicted  by  the  conventional  two-step  7,eldovich  mechanism  when  oxygen  atoms  are 
assumed  to  be  in  equilibrium  with  oxygen  molecules,  hut  are,  in  general,  less 
than  the  predicted  rates  when  oxygen  atoms  are  assumed  to  be  in  dynamic  steady 
state  with  the  CO-air  reactions. 

9.  Circumstantial  evidence  is  offered  in  support  of  the  proposition  that 
the  concentrations  of  nitrogen  oxides  in  the  CO-air  flames  are  near  their  equi¬ 
librium  values  for  pressures  above  about  20  atm.  The  case  for  equilibrium  is 
not  clear  at  the  lower  pressures,  but  it  is  suspected  that  kinetics  probably 
influence  or  even  control  the  formation  of  nitrogen  oxides  at  the  lower 
pressures. 
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APPENDIX  A 


METHOD  OF  CALCULATING  PSEUDOEQUIVALENCE  RATIO 

The  pseudoequivalence  ratio  <•>'  was  obtained  from  the  composition  of  the 
effluent  from  the  high-pressure  chamber.  It  was  calculated  from  its  definition 


^  =  <nCO/rVact  =  (WX 02)act 
(nco/no2)stoich  (xC0/x02)stoich 


where 


xco\ 

*°2/stoich 


The  ratio  of  CO  to  02  in  the  numerator  was  obtained  from  an  analysis  of  the 
effluents  as  the  ratio  of  the  unburned  CO  plus  the  burned  CO  (C02)  to  the  sum 
of  the  excess  oxvgen  remaining  after  combustion  and  the  amount  of  oxygen  con¬ 
sumed  during  the  combustion,  as  indicated  in  the  following  equation: 


Q 


act 


XC0  +  XC02 


x02  +  -  x002j 


eff 


The  air  flow  rate  was  calculated  from  <!>'. 
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APPENDIX  B 


THERMOCOUPLE  CORRECTIONS  FOR  RADIATION  LOSSES 

■to  correct  the  raw  thermocouple  readings  for  radiation  losses,  the  follow¬ 
ing  formula  of  Kasxan  from  reference  20  was  used: 


1 .25ead3/4/u  V  /4 

Tf  -TC-- ■—  {-)  (V-V, 

The  value  of  e  was  taken  to  be  0.25  (refs.  2i  and  22),  d  was  obtained  from 
micrometer  measurement  of  the  thermocouple  bead,  the  product  Pv  was  obtained 
from  the  expression  pv  =  m«j./A,  and  X  and  u  were  evaluated  at  Tf  for  a 
CX>2  and  N2  mixture  resulting  from  the  stoichiometric  combustion  of  CO  and  air. 
Mixture  properties  for  (J  and  X  were  calculated  using  Wilke's  equation 
(ref.  23)  and  Brokaw's  equation  (ref.  23),  respectively.  The  pure-gas  proper¬ 
ties  for  p  and  X  were  obtained  from  reference  24.  The  quantity  Tw4  was 
considered  negligible  compared  with  Tc4,  and  was  therefore  dropped  from  the 
equation. 
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ESTIMATION  OF  FLAME  TEMPERATURE  AS  A  FUNCTION  OF  PRESSURE 

An  approximate  estimate  of  flame  temperature  was  obtained  by  mathemati¬ 
cally  accounting  for  those  sources  of  heat  loss  from  the  flame  which  reduce  its 
temperature  to  some  value  below  the  adiabatic  flame  temperature,  as  follows. 

When  there  is  no  heat  loss  from  the  flame. 

Heat  Sensible  Sensible 

generated  +  heat  in  *  heat  out 

where  the  heat  generated  is  from  the  heat  of  combustion  and  the  sensible  heat 
refers  to  the  heat  content  of  the  gases  flowing  into  the  flame  zone  (CO  and 
air)  and  out  of  the  flame  zone  (CO2  and  N2) .  Stated  mathematically,  this 
becomes 


^RxmCO  +  ^  cp,  inni  (Tin  ”  Tref)  =  },  cp,ono(Tad  ~  Tref^ 


(Cl  ) 


When  there  is  heat  loss  from  the  flame. 

Heat  Sensible  Sensible  Heat  radiated  Heat  radiated  Heat  convected 

generated  +  heat  in  =  heat  out  +  to  chamber  +  to  burner  +  to  burner 

or 

^RxmCO  +  I  cp,  innin  ^Tin  ~  Tref'  =  5"  cp, o^^f  “  Tref)  +  ^ftfO(Tf4  ~  Tw4) 

+  Abrfd(Tf4  -  Tb4)  +  hAb(Tf  -  TV  (C2) 

Equating  equations  (C1 )  and  (C2),  collecting  terms,  neglecting  Tw4  and  Tb4 
relative  to  Tf4,  and  approximating  Tb  by  Tjn,  gives 

l  cp,0no(Tad  -  Tf)  =  (Af  +  Ab)efaTf4  ♦  hAb(Tf  -  Tin)  (C3) 

The  term  ^  cP(0n0  was  evaluated  for  a  CO2  and  N2  mixture  resulting  from  the 
stoichiometric  combustion  of  CO  and  air  and  from  the  known  flow  rate  for  CO. 

The  quantities  Af  and  Ab  were  obtained  from  measurements  of  the  external 
boundaries  of  the  luminous  zones  in  the  color  photographs  in  figure  4  and  are 
shown  plotted  in  figures  16  and  17.  The  emissivity  €f  was  obtained  using  the 
method  of  reference  25.  For  this  purpose,  the  luminous  zones  in  figure  4  were 
taken  to  consist  of  a  OO2  and  N2  mixture  resulting  from  the  stoichiometric  com¬ 
bustion  of  CO  and  air,  the  radiation  path  length  L  was  obtained  from  measure¬ 
ments  of  the  thickness  of  these  luminous  zones  in  the  upper  part  of  the  flames 
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(L  is  shown  plotted  in  fig.  18),  and  corrections  for  pressure  broadening 
of  the  0C>2  bands  were  made  by  employing  correction  factors  obtained  from 
reference  25. 

At  1  atm,  where  Tf  is  known  from  direct  thermocouple  measurement, 
equation  (03 >  was  solved  for  the  heat-transfer  coefficient  h,  giving 
0.362  cal/cm^-min-K  for  the  1 46-sccm  CO  flame,  and  0.331  cal /cm‘-mi n-K  for 
the  73-sccm  CO  flame.  Because  the  heat-transfer  coefficient  is  a  function 
of  the  mass  flow  rate,  the  various  transport  properties  of  the  gases,  and  the 
geometry  of  the  burner,  and  because  these  are  largely  invariant  with  pressure, 
these  calculated  heat-transfer  coefficients  are  taken  to  be  applicable  at  all 
pressures.  This  procedure,  ther,  allows  the  solution  of  equation  (03)  for  Tf 
at  all  other  pressures.  The  resulting  temperature  curves  are  shown  plotted  in 
figure  8.  Because  of  the  extensive  calculations  and  necessary  approximations 
involved  in  obtaining  these  temperature  curves,  they  must  be  considered  approx¬ 
imate;  nevertheless,  they  are  believed  to  represent  a  fairly  reasonable  trend 
of  flame  temperature  with  increasing  pressure. 
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APPENDIX  D 


ESTIMATION  OF  THEORETICAL  CO  REACTION  RATES 

Howard,  Williams,  and  Fine  (ref.  7)  give  for  the  global  rate  of  oxidation 
of  CO  to  form  CO2,  in  the  presence  of  water  vapor 


.G  -  1 . 3  «  1  O'  «  5098/T  Ccoci^cig  (01 1 

In  a  diffusion  flame,  combustion  within  the  flame  zone  may  reasonably  be 
assumed  to  occur  at  an  overall  average  equivalence  ratio  4>  =  1  ;  hence,  for 
the  present  CO-air  flames,  the  initial  composition  can  be  taken  to  be 


xCo  =  0.296 
xq2  =  0.148 
xN2  =  0.556 

Hydrogen  is  almost  always  present  in  small  concentrations  because  of  virtually 
unavoidable  contamination  with  water  vapor,  or  by  intentional  addition.  For 
most  of  the  experimental  runs,  the  CO  cylinder  with  99.3  ppm  H2  was  used  since 
no  difference  in  results  was  observed  between  this  concentration  and  the  lower 
concentration  of  10  ppm  H2  in  the  second  cylinder.  Therefore,  for  purposes  of 
estimating  the  theoretical  CO  combustion  rates,  the  concentration  in  the  fuel 
gas  was  taken  to  be  approximately  100  ppm  or  xH2  =  10"4.  In  a  stoichiometric 
mixture  of  CO  fuel  and  air,  this  gives 

xH2  =  2.96  x  10~5 

Taking  xH2q  -  xH2,  and  substituting  the  given  values  of  mole  fractions  into 

equation  (Dl )  at  T  *  Tf,  the  rate  curves  labeled  rG  in  figure  10(a)  are 
obtained. 

Several  other  expressions  for  the  global  rate  of  oxidation  of  CO  also 
appear  in  the  literature.  (See,  for  instance,  refs.  26  to  29.)  For  the  same 

1 

flame  conditions,  these  expressions  predict  rates  from  as  much  as  6—  times  the 

2 

rG  values  in  figure  10(a)  to  as  low  as  1/5  these  values,  and  some  have  somewhat 
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different  temperature  dependencies.  The  rate  expression  of  Howard,  Williams, 
and  Fine  (ref.  7)  was  selected  here  because  it  is  more  recent  and  because  it  is 
based  on  more  extensive  data.  In  addition,  it  predicts  rates  which  are  reason¬ 
able  averages  of  those  predicted  by  the  other  rate  expressions. 

An  alternate  estimate  of  the  rate  of  combustion  of  CO  can  be  obtained  from 
the  elementary  reaction  steps 

CO  +  OH  -*■  CO2  +  H  (d2) 

CO  +  O2  -*■  CO2  +  0  (d3) 

CO  +  0  +  M  -*■  C02  +  M  (D4) 

The  reaction  rate  constants  kf  associated  with  these  reactions  are  (ref.  8) 

k2  =  4  x  1 O1 2  e-4026/T 
k3  =  1 .6  x  1  o1 3  e-20634/T 
k4  =  6  x  IQ13 


with  consistent  units  in  moles,  cm3,  and  sec.  Reaction  rates  r2«  r3,  and  r4 

associated  with  reactions  (D2)  to  (D4)  were  evaluated  for  both  the  1 46-sccm  CO 
flame  and  the  73-sccm  CO  flame  as  follows.  The  hydroxyl  radical  and  oxygen-atom 
concentrations  for  use  in  reactions  (D2)  and  (D4)  were  determined  for  equilib¬ 
rium  conditions  at  T  =  Tf  and  for  an  equivalence  ratio  of  hydrogen  appropri¬ 
ate  to  the  present  system  of  100  ppm  H2  in  the  CO  fuel,  <(>H2  =  1  0-4  (i.e., 

when  the  fuel  is  taken  to  be  at  an  overall  average  $  of  unity  with  respect  to 
air).  The  results  were  then  scaled  for  xq2  =  0.148  in  the  flame.  The  rate 

r3  was  calculated  using  the  mole  fractions  for  CO  and  O2  given  earlier.  Com¬ 
parison  of  these  calculated  rates  showed  that  r4  <<  r3,  and  r2  is  roughly 
10  percent  of  r3  for  the  1  46-sccm  CO  flame,  and  roughly  2  percent  of  r3  for 
the  73-sccm  CO  flame.  Hence,  a  rough  estimate  of  the  theoretical  reaction 
rates  for  CO  can  be  written 

rT  =  r2  +  13 

where,  for  the  1  46-sccm  00  flame,  r>p  *  1.1 r3,  and  for  the  73-sccm  00  flame, 
rT  •  1.02r3.  These  rates  were  calculated  at  T  =  Tf  and  <}>  =  1  and  are 
plotted  in  figure  10(a).  The  close  agreement  between  r<p  and  rg  for  each 
of  the  two  CO  flow  rates  lends  confidence  that  these  rate  predictions  are  of 
the  proper  order  of  magnitude. 
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ESTIMATION  OF  DIFFUSION  RATES 

To  estimate  the  rate  of  diffusion  of  CO  from  the  core  of  the  flame  into 
the  reaction  zone,  Fick's  law  can  be  written  as 


^  dcCO 

J  =  -D  - 

dfi. 


where 

-► 

J  diffusional  flux  of  CO 


D  diffusion  coefficient  for  CO  (taken  as  CO  into  N2) 


rtr00 

-  concentration  qradient 

dH 

Expressinq  Cro  in  terms  of  x<-^>  by  the  perfect  qas  law,  and  wrttinq  the  con¬ 
centration  qradient  in  finite  difference  form,  there  results 


>  p  (xCO^core  “  *xCO^edqe 

J  =  D - 

RTf^  L 


where  the  subscripts  "core"  and  "edqe"  refer  to  the  inner  and  outer  edqes  of 
the  flame  zone,  respectively,  and  I,  is  the  thickness  of  the  flame  zone  (the 
same  as  the  radiation  path  lenqth  in  appendix  C) .  Takinq  (*00) core  ”  1  an<1 
(xCo)edqe  =  0,  and  writinq  the  diffusional  flux  as  a  volumetric  diffusional 
rate  rp  usinq  the  thickness  of  the  flame  zone  to  define  the  diffusional 
volume,  there  is  obtained 


rD 


PP 

RTfL2 


Values  for  D  were  calculated  usinq  the  method  of  Hirschfelder  as  described  in 
reference  23,  and  L  was  taken  from  fiqure  10. 
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ESTIMATION  OF  OXYGEN-ATOM  CONCENTRATION 

In  the  combustion  system  of  CO  and  air,  the  primary  reactions  involving 
oxygen  atoms  are  (see,  for  instance,  ref.  8) 


CO 

+  o2  -*•  co2 

+  0 

(Fl  ) 

CO 

+  0  + 

M  +  C02  +  M 

(F2) 

°2 

+  M  * 

20  + 

M 

(F3) 

0  -* 

•  n2  - 

NO  + 

N 

(F41 

N  -4 

■  02  - 

NO  + 

0 

(FS) 

The  reaction  rate  constants  k4,  in  units  of  moles,  cm3,  and  sec,  for  these 
reactions  are  given  by 

kT  =  1  .6  x  l  o1  3  e_20634/T 
k2  *  6  x  TO13 


from  reference  8,  and  by 

k3  =  3.56  x  1018  -  e-5938fi/T 
T 

k3r  -  I  x  TO14 

k4  =  1.36  v  1014  e-37947/T 

k5  =  6.4  x  TO9  Te-3M5/T 

from  reference  9,  where  all  rate  constants  are  for  the  forward  reaction  as 
written,  except  the  rate  constant  designated  with  a  superscript  r,  which  is 
for  the  reverse  reaction.  Since  C02  and  NO  are  stable  reaction  products  which 
exit  the  flame  zone,  only  the  forward  rates  of  reactions  (Fl),  (F2),  (F4), 
and  (F5)  are  considered;  but  since  oxygen  atoms  are  not  a  stahle  reaction 
product,  the  reverse  rate  of  reaction  (F3)  is  considered  in  addition  to  its 
forward  rate.  Employing  a  steady-state  assumption  for  oxygen  atoms  in  the 
flame,  the  rate  of  production  of  oxygen  atoms  is  set  equal  to  the  rate  of 
destruction  of  oxygen  atoms.  This  equality  can  be  expressed  in  terms  of  the 
rates  of  reactions  (Fl )  to  (F*>)  as 
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r?  +  2r 3  +  r«j  «  r2  +  2r3r  +  T4  (F6) 

where  all  rates  r  stand  for  the  forward  rates  of  the  reactions  as  written, 
except  the  rate  with  the  superscript  r,  which  designates  the  reverse  reaction. 
An  order  of  magnitude  analysis  of  equation  (F6)  for  temperatures  between  1600  K 
and  2200  K,  for  all  pressures  between  1  and  50  atm,  and  for  mole  fractions 
xCO  ”  0.296,  Xq2  ”  0.148,  xNz  «  0.556,  with  xM  -  1.000,  and  with  the  assump¬ 
tion  that  atomic  nitrogen  is  in  equilibrium  with  molecular  nitrogen,  indicates 
that 


H  »  2r  3 
ri  »  r5 

r2  *  0(2r3rl  for  high  oxygen-atom  concentrat ions 
r 2  >>  r 4 

Hence,  to  a  reasonable  approximation  at  all  temperatures  and  pressures  in  the 
flame, 


rl  “  r 2  +  2r3r 


<F7> 


Substitution  of  the  appropriate  rate  expressions  into  equation  (F7)  permits  a 
solution  for  the  steady-state  concentrat ion  of  oxygen  atoms.  Tn  terms  of  mole 
fractions,  this  concentration  is 


x0  =  -0.1  5x^0  +  ^0.0225x^0  6.565  -  ,,-20634/T  x^o^ 

Taking  x^  =  0.296  and  x0  =  0.148,  this  simplifies  to 


x0  =  -0.0444  +  [1.971  x  10~3  +  0.2876  -  e-20634/T 

V  p 


J/2 


(F8) 


(F9) 


where  the  subscript  ss  is  appended  to  designate  that  this  oxygen-atom  concen¬ 
tration  in  not  an  experimentally  determined  value,  but  is  instead  an  estimated 
value  derived  from  the  present  steady-state  analysts.  It  should  he  pointed  out 
that  since  the  reverse  of  reaction  (Fi )  was  neglected  in  the  analysis,  the  val¬ 
ues  of  x0  rr  obtained  from  equation  (F9)  probably  represent  an  upper  limit 
for  Xq  which  might  be  expected  in  the  flame.  The  (F2)  reverse  reaction, 
which  produces  oxygen  atoms,  was  also  neglected?  hut  because  this  reaction  is 
highly  improbable,  its  contribution  to  the  concentration  of  oxygen  atoms  is 
insignificant.  Values  of  x0> evaluated  from  equation  (F9)  at  the  appropri¬ 
ate  values  of  Tf  and  p,  are  plotted  in  figure  12  as  x0>  ss/^oleq^* 
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Figure  2.-  Schematic  diagram  of  burner  and  sample  collection  system. 
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Figure  4.-  Photographs  of  flame. 
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